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Abstract 

We present a propagator formalism to investigate the scattering of photons by 
a cavity QED system that consists of a single two-level atom dressed by a leaky 
optical cavity field. We establish a diagrammatic method to construct the 
propagator analytically. This allows us to determine the quantum state of the 
scattered photons for an arbitrary incident photon packet. As an application, 
we explicitly solve the problem of a single-photon packet scattered by an 
initially excited atom. 
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I. INTRODUCTION 



The engineering of novel quantum states of photons is a topic of interest fundamentally 
and for applications. As is well known in atomic and particle physics, exotic states could 
result from scattering processes. This suggests light scattering by suitable quantum systems 
can be an important tool for state engineering. A common and yet realistic situation in 
quantum optics involves the scattering of photon wave packets, comprising a few photons, 
from atoms and other photons situated inside an optical cavity. The atoms and the cavity 
field together form a cavity QED system that scatters incident photon packets (Fig. 1). The 
main question is: What is the quantum state of the scattered photons? As the cornerstone of 
a series of investigations of this topic, we establish in the present paper a propagator method 
to describe the scattering of photon packets from a cavity where a two-level atom, either 
in its ground or excited state, is placed with the companion of some quasi-mode photons. 
After the interaction, the input photons and the quasi-mode photons as well leak out of the 
cavity and become the output state of the system. The propagator method proposed here is 
capable to extract the details of the output state and thus form a powerful tool to analyze 
the physical content underlying the scattering process. 

In addition to its potential applications in the study of novel quantum states, the model 
considered here is in fact a generalized variant of the Jaynes-Cummings model (JCM) [|T|, 
which has been a fundamental model in the realm of photon-atom interaction ever since its 
first introduction 0. Not only does it provide the simplest description of an atom interact- 
ing with quantized fields in optical cavities, JCM is also an important tool for controlling 
quantum states . One prominent generalization of the JCM is the inclusion of the leakage 
effect of the cavity, which is inevitable in any real experimental setup. Indeed recent studies 
in cavity QED have also emphasized the important role of dissipation effects in quantum 
information devices [Q. 

In conventional approaches to JCM in a leaky cavity it is customary to consider the 
space outside the cavity — the environment — as a Markovian bath, followed by solving 
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the master equation of the reduced density matrix of the cavity field Thus, energy 

inside the cavity flows unidirectionally to the environment, which plays the role of a sink of 
energy and the information on the exact state of the environment is reckoned unimportant 
and is lost during the evolution. 

The model discussed in this paper distinguishes itself from other dissipative JCM's in 
that the cavity and its environment play an equal role during the interaction. These two 
parts, constituting the "universe" in consideration, communicate their information through 
partially transmitting walls dividing them. Photons inside the cavity, described in terms 
of discrete quasi-modes, interact with the atom and eventually leak out of the cavity. On 
the other hand, incident photons characterized by continuous wave number enter the cavity, 
participate in the interactions with the atom, and leave the cavity in its flnal state as well. 
More importantly, these two kinds of photons lose their respective identities after entering 
(or returning to) the environment and interfere to form novel photon states. It is for this 
reason that we devote the current paper to a proper formalism describing this generalized 
JCM model. 

To this end, we present here a propagator method that solves directly the quantum state 
of the whole amalgamated system, consisting of the atom, the cavity and the environment. 
Rather than discrete cavity modes, the atom now couples with the continuous fleld modes 



of the whole system ||TO|-[T^. In this sense, the artiflcial separation between the cavity and 
its environment is completely eliminated, and the effect of dissipation can be treated in 
a more rigorous and fundamental way. Once the propagator of the system is known, the 
exact details of the output photons can be obtained immediately without further ado. This 
knowledge is bound to be important for issues related to quantum information and quantum 
measurement. 



The pure state approach has previously been used to study resonance fluorescence in 



and spontaneous atomic decay in a vacuum cavity, and exact solutions have been found 
|jn|,[TB|-|T7| . The major objective in the present paper is to formulate a comprehensive theory 



studying the dynamical response of a two-level atom under the influence of photons inside 



the cavity and from the environment. In a recent paper by the authors |T^, a less general 
case with the input photons restricted to the cavity mode (i.e., quasi-mode) was considered 
and solved by method of Laplace transform. This restriction is now removed in this paper. 
By evaluating the propagator of the system, whose construction is based on a diagrammatic 
expansion method, one can readily handle situations with an arbitrary initial state. Under 
the rotating-wave-approximation, the full propagator is a block diagonal matrix with a se- 
quence of 2 X 2 matrices forming its diagonal block. Each of these matrices is a propagator of 
the system in a subspace categorized by the total excitation number (to be defined rigor- 
ously in Sect. IV), providing a natural classification of relevant propagators. A construction 
method is established here to express the A^-excitation propagators in terms of those of the 
lower excitations. Thus, the full propagator of the system can be obtained in a systematic 
way, paving the way for further investigation into the quantum state of the output photons. 

This paper is organized as follows: In Sec. II, we describe the system under study, its 
normal modes, the atom-field interactions and obtain the Hamiltonian. An introduction on 
the propagator method and the related diagrammatic expansion will be given in Sec. III. 
The general solution for the propagator is derived in Sec. IV. We then apply the propagator 
for A^ = 2 to study the scattering of a single-photon packet from a two-level atom in its 
excited state in Sec. V. We draw our conclusion in Sec. VI. 



II. THE MODEL 
A. The cavity and continuous modes 

We consider here a two-sided Fabry-Perot cavity in one-dimensional space. The two 
transmitting mirrors (hereafter referred to as mirrors R and L) are placed at x = ±L/2. 
The set of continuous modes of such a Fabry-Perot cavity is known in the literature |[T0| , p!3| . 
In the following, we will briefly sketch the main result for reference. 

The two mirrors enclosing the cavity are modeled by two thin dielectric slabs of thickness 
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/ and refractive indices {a =L,R). In the limit / and ^ oo such that n^/ 

is finite, the corresponding complex amplitude refiection and transmission coefficients ta 

(a = L, R) are in turn expressed as 



tn 



2 7/k jjify 
2 



(2.1) 
(2.2) 



2 %}\,^(y^ 

where k is the wave number. 

There are two independent sets of modes for the entire system (polarization is ignored in 
the current study), namely, the left-propagating modes Ui{k,x) and the right-propagating 
modes Ur,(A;,x), which are given by |10|,|13[: 

e'^^ + Ri,{k)e-'^'' -oo < x < -L/2 

h{k)e''''' + JL{k)e-'^'' -L/2 <x < L/2 
TL(A;)e^'=^ L/2 < x < oo 



UL{k,x) 



(2.3) 



UR{k,x) 



TR(A;)e-^*^^ 
jR(fc)e-^^^ + jR(A;)e^'=^ 
e-^'^^ + i?R(fc)e*'=^ 



-oo < X < -L/2 
-L/2 <x <L/2 
L/2 < X < oo 



(2.4) 



where 



R^(k) = {rLC"^'^^ + rRe^'=^+2*'^''s*L}/D(A;) 
h{k) = tjD{k), 
Mk) = t^rne''''^/D{k), 
n{k) = Tn{k)=t^tn/D{k), 



(2.5) 
(2.6) 
(2.7) 
(2.8) 



with 



D{k) = 1 - TLrRC 



2ikL 



(2.9) 



Likewise, /r(/i;), Jr(A;) and RR{k) can be obtained by interchanging the roles of L and R in 
the above equations. The left-luminating modes are shown in Fig. ^ The mode functions 



ui,{k,x) and u^_{k,x) {k > 0) together form a complete orthonormal set in — oo < x < oo, 
satisfying the orthonormal condition: 

/oo 

n{x)ua{k, x)u*p{k', x)dx — 2'K5ap5{k — k') , (2-10) 
-oo 

where a, (5 =L,R. 

The quasi-mode frequencies, km {m — 1, 2, 3, . . .), defined by the zeros of D{k), are given 
explicitly by 

T ^ 1 / \ ?717r ,^ ,,^ 

km = ^HrLTR) + — . (2.11) 

For real-valued k close to a quasi-mode frequency, i.e., |A; — km\ <^ tt/L, we have 

D{k) = -2iL{Ak + m^) , (2.12) 

where Ak = k — kc, with 

7T17T 1 

kc = —- ^Im[ln(rLrR)] , (2.13) 

«c = -^Re[ln(rLrR)] . (2.14) 

Here Hc is the decay rate of the cavity. 

B. The two- level atom and interaction Hamiltonian 

Consider a system of a two-level atom placed inside the cavity described above (or other 
leaky cavities in general) at x = xq- The ground-state energy of the atom is arbitrarily taken 
as zero, while the excited-state energy is cua in units of ^ = c = 1. The full Hamiltonian of 
the system in the rotating-wave-approximation is given by 

H = ujAcr+cr- + / k{alj^akL + alJ^akR)dk 
Jo 

+ / {[9L{k)akL + gR{k)akR]a+ + h.c.}dk, (2.15) 
Jo 

where Uka and a^^, {a — L, R) are the annihilation and creation operators of the ka-mode 
photon, a± are the pseudo-spin flip operators of the atom. They satisfy the usual commuta- 



tion relations: 



(^ka, O'k'a' 



S{k — k')Saa' , [dka, O'k'a'] — ^L' ^k'a' — i^ka, <^± 



0, and {(7_,(T+} = 1. The coupling constant ga{k) of the atom with the ka-mode photon 
depends on the dipole moment of the atom and is proportional to Ua{k,XQ). Here we are 
particularly interested in the case where the transition frequency of the atom is close to 
one of the resonance frequencies of the cavity, say, kc- Hence, only those continuous modes 
with frequencies near kc have significant interactions with the atom and in this single-mode 
approximation, ga{k) oc {Ak + iKc)~^ ■ 

Despite that the atom ostensibly couples with both the R and L modes, it is always 
possible to use a unitary transformation to redefine the photon modes so that the atom 
interacts with only one set of modes. For concreteness, we consider in the present pa- 
per a symmetric cavity with identical mirrors and the atom being situated at its center. 
However, it can be proved that generalizations to cases with dissimilar mirrors, arbitrary 
atomic position and one-sided cavity are straightforward. For this specific model, we have 
Uhik, xq) — u^,{k, xo), ghik) — gnik). Accordingly, one can define a new basis of photons by 
the unitary transformation 

0'k+ = ;^("feL + akR) , (2.16) 

Ofe- = ;^("feL - OfeR) ■ (2.17) 

It is readily observed that the "— " modes do not couple to the atom. In the following 
discussion, we will ignore the "— " modes and will focus on the evolution of the modes, 
with the index being suppressed. Hence, the full Hamiltonian reduces to 

H = uuaO'+o'- + / kaiakdk + / [g{k)ak(J^ + g*{k)aia_]dk , (2-18) 
Jo Jo 

where ak here actually denotes ak+ and g{k) = \/2gR{k). 

III. PROPAGATOR AND FEYNMAN DIAGRAMS 

The dynamical response of a system with a Hamiltonian H is governed by the retarded 
Green's function that satisfies the equation: 
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(^^±-H^K^(t,t') = 6{t-t'), (3.1) 

and is null for t < t' . Obviously, Kj^ can be given explicitly by K+{t,t') = e-'"^^-^'^9{t -t'), 
where 6{x) is the Heaviside step function. Accordingly, its Fourier transform, G{uj), defined 
by the relation 

/oo 
K+{t,0)e''''dt (3.2) 
-oo 

and termed the retarded propagator G{uj), can be expressed symbolically as 

Giu) = -1- , (3.3) 

LU — ri 

where the prescription u = lim (u + ie) is assumed hereafter. 

To establish a diagrammatic expansion for G{uj), we separate the Hamiltonian in the 
form: H = Hq + V, where 

Hq = u!A<y+<y- + / kaiakdk (3.4) 

JO 

is the free atom-field Hamiltonian, and 

V= / [g{k)ak(x+ + g*{k)ala.]dk (3.5) 





represents the interaction between the atom and the field. Hence, the propagator can for- 
mally be expanded in a power series of V, yielding 



LO — Hq — V LJ — Hq UJ — Hq UO — Hq 

+^7r^^7r^^7r + (^.e) 

U — Hq UO — Hq UO — Hq 

The propagator so defined is an operator. In the energy-eigenstate basis, the transition 
amplitudes are given by the matrix elements of the propagator, also referred to as "prop- 
agator" hereafter. The main purpose of this paper is to calculate these amplitudes by the 
associated Feynman diagrams as illustrated in Fig. |^. The basic construction rules and 
interpretations of these diagrams are specified as follows: 
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1. External and Internal Lines: Photons are represented by wavy lines labelled by their 
momenta k. Atoms in the excited and the ground states are respectively represented by 
solid and dashed lines . 

2. Vertex Factors: Each vertex contributes a factor g{k) (photon absorption) or g*{k) 
(photon emission) to the associated amplitude. 

3. Free Propagators: Each segment in between successive vertices contributes a factor 
{uj — E)^^ to the associated amplitude, where E is the energy of the free Hamiltonian (i.e., 
Ho) in this segment. 

4. Integrate Over Internal Momenta. For each internal momentum k, write down a factor 
dk and integrate. 



IV. EVALUATION OF THE PROPAGATORS 

In this section, we evaluate the matrix elements of the propagator using the diagrammatic 
rules stated in Sec. Owing to the rotating-wave-approximation, it is readily observed 
from Eq. ( p.l8| ) that the total excitation number 

roo 

N = a+a_ + / alakdk (4.1) 
Jo 

is a constant of motion, resulting in vanishing propagators from an initial state to a final 
state with different excitation numbers. Therefore, the full propagator can be represented 
by an infinite sequence of 2 x 2 matrices, each characterized by its excitation number A^. 
Hereafter we will, for convenience, define an energy eigenstate of Hq with an excitation 
number N by: 

\p; ki, k2,--- kN-p) = , ^ ,, 4i42 ■ ■ • «L_p^+|0; 0) = b; K^-p) ) (4-2) 

^{N-p)\ 

where p = 0(1) if the atom is in its ground (excited) state, and |0) is the vacuum-field state. 
The factor 1 / ^ {N — p)\ is introduced here to take care of the multiple-count of the bosonic 
states in integrations and \p; ki, k2, ■ ■ ■ k^-p) is not necessarily normalized to unity. 
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In terms of this notation the four propagators with excitation number N are 

Gl^\u;; K^_p, K'^-,) = (p; K^-p| \q; K'^-,) , (4.3) 

U) — n 

where p, g = 0, 1. As an example, the propagator of zero excitation number, governing the 
propagation of a ground-state atom in vacuum, is trivially given by 

GS(a;;0,0) = -, (4.4) 

UJ 

for there is only one diagram, namely, the free propagation of the collective vacuum state. 



A. Quasi-mode propagators 

Before proceeding to explicit evaluation of general propagators, we introduce here the 
concepts of quasi-mode photon states and quasi-mode propagators. A normalized quasi- 
mode single-photon state is defined by 

lie) =410), (4.5) 

where a\ is the effective creation operator for the quasi-mode: 

al^^j^^ dkg*{k)a\. (4.6) 

The quantity A in the normalization constant is the coupling strength defined by 

A=/ g{k)g*{k)dk. (4.7) 
Jo 

Physically speaking, the state |lc) is, in a perturbative sense and also in the single- mode- 
approximation, the cavity field set up by the atom during de-excitation process. Similarly, 
the atom-field state with excitations, where there are N — p quasi-mode photons and p 
atomic excitation, is defined by 

\P-AN-P).) ^ , ^ (at)(^-^)aj^|O;0). (4.8) 
'{N -p)] 



We therefore accordingly define the A"-excitation quasi-mode propagator by 
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$(f M = {p- {N - p),\ -— ^ \q; {N - g),) . (4.9) 



From the Hamiltonian given in Eq. ( p.l8|) , the importance of the quasi-mode propagator is 



readily clear. Aside from the input and output photons, those present in the intermediate 
states are all quasi-mode photons. Thus, these quasi-mode propagators form the backbone 
of our theory from which other propagators can be derived. 

From definition ( [4.9| ), it is obvious that the vacuum propagator in Eq. ( |4.4| ) is the simplest 



quasi-mode propagator: ^'qq{uj) = Gqq {u; (j),(f)) = u ^. We begin with the quasi-mode 
propagator of single excitation 

$SVH = (1;0|^--|1;0). (4.10) 

00 — 12 

All relevant diagrams are shown in Fig. ^ and the physical picture indicated by the diagrams 
is clear. The excited-state atom may freely propagate in vacuum, possibly followed by equal 
numbers of emissions and absorptions of quasi-mode photons, and exits in its excited state. 
In Fig. H, bold-wavy lines are used to represent intermediate quasi-mode photon state, to 
distinguish it from the input and output normal-mode states. 

According to the Feynman rules, ^^i{uj) is given by an infinite series: 

^^i^H = ^—(i + ^^+ . + , + , (4.11) 

CO — loa \ CO — loa [io — uja) — ooa) J 

where 

cH=r~£M,,. (4.12) 

JO 00 — K 

Adopting the single-mode approximation and, as usual, extending the lower limit of all the 
/c- integrations from to — oo, we find 

CH = 7^-—. (4.13) 

oo — kc + IKc 



Therefore, Eq. ( 4.11 ) can be expressed in the closed form 

oo-n^l'^^ Lo- n^^^ 



^lli^) = ^-m + ^Jl^^ (4-14) 
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where 



1 ± 



^ = ^ + ( iV 



{ua - kc + mc)/2 
{ujA -K + mc)V4 + NX 

1 



(^c if^c) i 



/ oja ~ kc + 
2 



+ NX 



(4.15) 



(4.16) 



for = 1, 2, 3, . . .. It is then obvious that V NX essentially plays the role of the A^-photon 
Rabi frequency. 

We have derived in the Appendix the quasi-mode propagator for an arbitrary excitation 
number A^. They are given by: 



A 



(N) 



+ 



A 



(N) 



1 - ) 1 - aL^) 



Vnx 



io - uja- {N - l){kc - iKc) 

{N) 



(4.17) 

(4.18) 

(4.19) 
(4.20) 



B. Propagators of single excitation {N = 1) 

With the help of the quasi-mode propagators, we can derive the simplest propagators 
for the case A^ = 1 to manifest the techniques in calculations aided by the diagrammatic 
scheme. We begin with the propagator G^ii{u!; 0, 0) = (1; (f)\{uj — H)^^\l; 0), which coincides 
with the first-order quasi-mode propagator given in Eq. (|4.14| ), i.e., 



Gff(o;;0,0) = <l>SlHo;). (4.21) 

One can, of course, follow the same route in the deviation of G^ii{uj; 0, 0) to obtain the 
other three propagators of single excitation. However, it is clear that once a particular A^- 
excitation propagator is derived, other members in the same class can be easily obtained in 
an alternative way by relating them to the one already obtained, as shown in Fig. ^. For 
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example, from Fig. ^ (a), the propag ator GS'{uj;k,(f)) = {0; k\{iu - H)~^\l ; 6) can be related 
to G'Sl^(a;;0,0)by: 

cSico; k, 0) = ^gS1)(u;; 0, 0) . (4.22) 

Likewise, we can show from Fig. ^ (b) that the propagator 

G^Siu;; 0, k) = (1; 01^— :|0; k) = ^Gff (.;; 0, 0) . (4.23) 
00 — H UJ — k 

Eqs. ( |4.22| ) and ( [4.23| ) also reveal a useful relation of the propagators, namely, the propagator 
G^^\caj; K^v-p, KV-g) can be obtained from G^J^\uj; KV-g, K^v-p) by simply replacing each 
g{k) with g*{k), and vice versa, which is obvious from the diagrammatic scheme and the 
Feynman rules. 

Finally, for the scattering of a photon from the k'-th mode to the fc-th mode by the 
ground-state atom, the corresponding propag ator is A;,A;') = {0;k\{iu - Hy^\0;k'). 

Similarly, from Fig. § (c), we have 

Qo k,k) = — + — T7^G\{ 0, . 4.24 

u — k' [ijj — k)[ijj — k') 



C. Construction of propagators of general excitation number N 

In general, an A^-excitation propagator can be calculated from those of lower excitation 
numbers. We will develop here a systematic approach to evaluate the propagators with 
N>2. 

First of all, in some diagrams there may exist "spectator photons" that do not inter- 
act with the atom in the whole evolution. These diagrams are said to be factorizable (or 
unlinked), and can be straightforwardly related to propagators of lower excitation numbers 
in a way to be stated explicitly in the following discussion. We therefore focus mainly on 
the linked diagrams that all the input and output photons take part in the interactions. 
These absorptions and emissions can take place in any order, resulting in different physical 
processes. 
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For each linked diagram we label it with a set of momenta, S = {ki, k2, ■ ■ ■ , k2N-p-q}, 
comprising the time-ordered momenta of the photons created or annihilated in the process. 
Here p and q are integers defined in Eq. (|4.3|) . The elements of S are taken from N—p output 
photons in the set Kat-p and N — q input photons in the set KV-g, with no repetitions. 
In other words, a particular S represents a particular sequence of the N — q absorptions of 
the input photons and N — p emissions of the output photons. For a given S, we introduce 
a quantity C^\uj]S) to denote the corresponding contributions to the propagator. The 
symbol "£" is used here to refer to "linked" diagrams. We note here that this convention 
should be accompanied by a final symmetrization of the propagator with respect to K.n-p 
and K'jv-g- 

There are at most only three kinds of photons present in any segment of the Feynman 
diagrams, namely, the input and output photons in continuous modes, and the quasi-mode 
photons. Between the 2N — p — q vertices associated with the sequence defined by S, 
the atom only interacts with the quasi-mode photons. Therefore the evolution is governed 
by the quasi-mode propagators derived. This explains the importance of the quasi-mode 
propagators discussed in Sec. ( [IV A| ). 

It is found that C^p^^{uj] S) can be written in the following compact form 

/N-q \ /N-p \ /2N-p-q \ 

4f ^) = ( n 3{k'i) \ I n 9\h) \ I n "^-si^ - 1 • (4-25) 

The physical picture of the above equation is clear. N — p times of absorptions of input 
photons and N — q times of emissions of output photons must take place at some time, with 
amplitude given by the first part of the equation. These events split the whole process into 
2N —p—q+1 segments. In a particular segment, the atom interacts only with the quasi-mode 
photons. This means that the propagator of this segment is governed by the quasi-mode 
propagators derived in the last section. The presence of continuous mode spectator photons 
in the i-th segment has the sole effect of shifting the frequency uj hy Ei, which is the total 
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energy of the continuous-mode photons present in the i-th segment Q. 

While not exphcitly shown, the C^\uj; S) in Eq. ([4.25| ) depends on S through the terms 
Ni, Pi, qi and Ei. 

Summing all diagrams corresponding to all possible S, we have 



(4.26) 



which is the propagator that includes all linked diagrams only. 

To include the unlinked diagrams, we go back to the case that some of the input and 
output photons act as spectators in the whole process. There may be one spectator, which, 
without loss of generality, can be assumed to be kiy_p (fcV-g)- The propagation of the 
remaining system is governed by A^^^^\ Similarly we can have two photons acting as 
spectators, which are assumed to be kjy^p ik' N~q) and k^-p-i (/cV-g-i), and so on. The 
maximum possible number of spectator photons is M = min{A^ — p, N — q}. Hence we have 

M 



G(f(o.;K^_„KV,) = EE 

sym j=0 



Y[ 5{kN-p-l - k'N-q-l) 
J=0 

i-1 



^^^^ ~ X! ^N-p-l] Kjv_j_p, rs. N-j-q) 
l=Q 



(4.27) 



In the above equation, denotes the symmetrization of the expression with respect to the 

sym 



input and output photons, and A is given by Eq. ([4.26|). 



^For example, in Fig. |6|, consider the first stage of the second term in the summation, the excited 
atom interacts with quasi-mode photon while the input normal-mode photon k'l acts as spectator 
photon. The evolution of the atom plus quasi-mode would be governed by the the single excitation 
quasi-mode propagator <I>^^^ {uj) in the absence of k'^. However when k'^ is present, the propagtor 
of the whole system in this segment becomes ( u; — k'^). 
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D. Example: Propagators with excitation number N = 2 

The general expression of the propagators is simple if the excitation number is not 
large. In this section we will derive the four propagators of excitation number N = 2. 

Consider first the propagator {u; ki, k[) = (l;A;i|(co' — H)^^\l; k'l) . According to 



Eq. ( 4.27 ), we have 

Gf^iu; k,, k\) = Kf}{u- ki, k\) + 5{kr - k\)K^S{^ - h'AA) , (4-28) 

where 

a!?(u;;A;i,A;'i) = ^4?(cu;S), (4.29) 
s 

with S = {ki, ki'} or {ki , ki}. One can readily show that (see Fig. 

dSiuj; S = {k,, k[}) = gik\)g*{k,)^^^^iuj - k,)<^?^ iuj)^^ (uj - k\) (4.30) 

and 

C^^^iu; S = {k[, k^}) = g{k\)g*{k,)^^h^ ' h) *S(^ - k, - (a; - k\) . (4.31) 

The total propagator is hence 

Gf^iu; ki, k\) = $ff (^ - fci) ( 5{k^ - fc/) + g{k\)g\k{)^^S{^ " ^'i) 

X 



' -7+ .. . . , .J. . , 



uj — ki — ki {oj — ki — kc + ii^c){'^ — W — kc + ii^c) 



(4.32) 



The other three propagators of the same excitation number can be evaluated simi- 
larly. However, as mentioned in Sec. |IV B| , these propagators can be obtained imme- 



diately from their relations with Gfi . For example, the propagator GQi{uj]kik2,k'i) = 
(0; A;i, k2\{uj - H)-^\l- k\) is given by (see Fig. (a)) 

^[.^(a;; k.k^. k\) = ^-—[g\k,)Gf^{uj- k^, k\) + g*{k2)Gf^{uj- k,, k\)] . (4.33) 

UJ — Ki — K2 
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Likewise, we can obtain the remaining two propagators with N = 2 from Figs. ^ (b) and 
(c), and the results are stated below: 

Gf^iu; h, k\k'2) = (1; A;i|^-|0; k\, k'2) 

to — n 

^ rWi)Gfl{u;- k,, k'2) + g{k'2)Gf}{uj- k,, k\)] . (4.34) 



LU — k'l — k'2 

G'S(^; kik2, k\k'2) = (0; fci, k2\^—\0; k\, k'2) 

uj — H 



^ -[5ik\ - k{)5{k'2 - k2) + 5{k\ - k2)5ik'2 - fci)] 



Lj — k'l — k'2 

+ Jl-^[9ik\)G^oii^; hk2. k'2) + gik'2)G^oii^; ^1^2, k'^)] . 

UJ ~ k'l — k 2 



(4.35) 



V. APPLICATION: SINGLE- ATOM SINGLE-PHOTON SCATTERING 

As an application of our method, we study the scattering of a single-photon wave packet 
by an excited two-level atom inside the cavity. In particular, we investigate how the spectral 
width of the incident photon affects the outcome of stimulated emission process. To begin, 
we consider the incident photon initially prepared in the "+" modes as defined in Eq. ( [^.16[ ). 
Hence the initial state is given by 



\^{t = 0)) = j dk'G{k')al,\l; 0) , (5.1) 

where G{k') is the photon amplitude. At a later time t, the state becomes 

m)) = J dkB{k, t)4|l; (j)) + l=JJdki dk2G{ku k2, t)4^4jO; 0) , (5.2) 

where the two-photon amplitude C{ki, k2) = C{k2, ki) satisfies the normalization condition, 

JJ \C{ki,k2)fdkidk2 = l. (5.3) 

For simplicity, we will only consider the resonance case: kc = wa- 

The long time state {ipit — > 00)) is determined by the asymptotic behavior of G{ki, k2, t). 
Utilizing the propagator obtained in Eq. (|4.33| ), we have 
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limC(ki,k2,t) = lim / due'''^' f dk'cf^ (u; kik2, k')C{k') 
The explicit form is given by 



(5.4) 



C{ki,k2,t-^ oo) -j=t 



-i{ki+k2)t I * 



9*ik^)^\^iki)Cik2) + g*ik,)g*ik2] 



X $SV(A;i)/i(A;i,A;2) + 
+{ki ^ /ca} , 



A 



— $ff(A;i)$ii^(A;i + A;2)/2(A;i,A;2 



.(2) 



(5.5) 



where 



h{ki, k2) 
h{ki, k2) 



hm 



- g{k')^^^^{ki + k2-k') 



ki-k' + iS 
9{k') 



C{k')dk' , 



<^\^{ki + k2-k')C{k')dk' , 



(5.6) 
(5.7) 



ki + k2 — k' — kc + iKc 

and {ki ^2} denotes the previous expression with ki and ^2 interchanged. It should be 
noted that the contour of the /i-integration should be closed in the lower half-plane. 

The main advantage of our formalism is that it can handle any state of the incident 
photon. As an example, we consider that C{k') is a lorentzian with a peak frequency equal 
to the cavity resonance frequency, i.e., 

1 



C{k') 



(5.8) 



Here the width is the spectral width of the incident photon. The pole is located in the 
lower-half-plane in order to ensure that the atom can only "feel" the photon for t > 0. 

In Fig. § we show the contour-plot of |C(A;i,A;2)P for A = O.Ik^, with various spectral 
widths of incident photons. When ^ 7sp (Fig. | (a)), the input photon has a very short 
pulse duration compared with the decay time of the atom. Therefore the incident photon 
is incapable of having sufficient interactions with the atom. As a result, the output state 
is approximately a direct product state of the input photon and the spontaneously-decayed 
photon of the atom, which corresponds to a "cross" shape in |C(A;i,fc2)p. A similar effect 
can be seen if the incident photon has a narrow width such that Kin ^ 7sp (Fig- 1§ (d))- In 
this case, the input photon has a very long duration and so it participates in the interaction 



mainly after the atom has reached the ground state, causing no interference with the photon 
emitted from spontaneous decay. 

However, interesting features show up when /tin is the same order as 7sp (Figs. § (b) and 
(c)). We see that the final two-photon amplitude is drastically different. For example, in 
Figs. ^ (b), there is an unexpected dip at the center and peaks at Ak ^ 7sp. In other words, 
although the frequencies of the input photon and the photon emitted in spontaneous decay 
both peak at kc-, it is very unlikely to have two photons with frequencies around kc in the 
output. Instead, the peak frequency has been shifted to Aki,Ak2 ~ \/nc- 

To understand the interference effects shown in Figs. ^ (b) and (c), we identify the 
contributions of relevant diagrams associated with the propagator. Fig. ^ (a) shows the con- 
tributions solely from the unlinked diagram in which the incident photon does not participate 
in the interactions at all. Obviously, the corresponding two-photon amplitude disagrees with 
the exact one. However, a much better agreement can be achieved if we include just the 
lowest order linked diagram given in Fig. ^ (b). The interference between the linked diagram 
and unlinked diagram produces a two-photon amplitude that is almost the same as the exact 
one (Fig. 11(c)). 

VI. CONCLUSION 

In this paper, we have developed a diagrammatic scheme to construct the propagator 
governing the interaction between an atom and photons in a leaky cavity. Under the as- 
sumption that the frequency dependence of atom-field coupling is a lorentzian, we found 
that the perturbation series is summable and so exact solution can be obtained for any 
excitation number N. The propagator provides an analytical tool to investigate the cavity 
QED effect for the photon scattering problem. In particular, the quantum state of output 
photons in continuous modes can be determined explicitly. 

Our results are illustrated by the derivation of propagators with excitation number = 
2, which are then applied to study the scattering of a single-photon wave packet from an 
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excited atom. We found that the spectral width of the incident photon can significantly 
modify the final two-photon amplitudes. This occurs when matches the cavity modified 
atom decay rate. Our calculations show that two output photons are entangled in the sense 
that their spectrum displays nontrivial correlation. With the aid of Feynman diagrams, 
we have identified the essential processes causing the interference. However, more detailed 
investigations are needed for a thorough understanding of the rich features in cavity QED 
assisted photon scattering problems. 
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CUHK4282/00P) and a direct grant (Project ID: 2060150) from the Chinese University of 
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APPENDIX A: DETAILS OF THE DERIVATION OF QUASI-MODE 

PROPAGATORS 

The main idea of the evaluation of the propagators is to move all the annihilation and 
creation operators to the right and left side of the perturbation series, employing the com- 
mutation relations, followed by noticing the fact that ak\(j)) — {(plal — and a"^ — — 0. 

The commutation relations between Uc, al with Hq and V are: 



ACKNOWLEDGMENTS 



al, V — —y/\a^ 

[Oc, V] = V^(7_ , 



+ ) 



(A2) 



(Al) 




(A3) 



(A4) 
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a. $ff (^) 

Consider the quasi-mode propagator 

$ffH = (0;iV,|^— |0;iVe). (A5) 

UJ — II 

The atom is initially at its ground state with photons in quasi-mode. The propagator 
measures the probability that the system remains in the same state finally. 

The interaction-free part can be evaluated by considering Eq. (|A3|) , which gives 

1 1 ^ r r 



1 



(A6) 



to — N(kc — i/tc) 

Paths with an odd number of interactions have zero contributions. And from the com- 
mutation relations, it is readily shown that 

V 10; Nc) = -^-^ TT / dkj g*(k.) at cr+|0; 0) Y , 



uj-Ho' ' ^' X^/^VW.f^J Y^-^c + mc-E,y^fc, 

(A7) 



and similarly 



^^;^)'|0;A^e) = yiJ^Y. {llj dk,g*{k,)al}j J dk g*{k)al\0-A) 



A 

X 



^ - Y^j^i kj - kc + iKc){uj - Y^j^i kj - uja) 
NX 



N 



n J dkjg*{kj)aljO; 



^ {uj-kc + - Ej=/ kj){uj -uj^- J2f=i^ kj) ' ''^^^ 



Hence 



(0; N,\^—V^—V^—\0- N,) 

UJ — H UJ — H UJ — H 

NX 

[uj — N{kc — inc)] [uj — uja — {N — l){kc — inc 



(A9) 
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It can be shown that the contributions of different paths with even times of interactions, 

1 f V V 

-1 i/2 



(0;iV, 



LO — Hq \iJ — Hq 
1 



Hence 



[a; - A^(A;c - ^c)] 



i/2+l 



NX 



uj - uja- {N - l){kc - iKc) 



(AlO) 



1 



tu — N{kc — inc) 
1 



where 



c 



[io - N{k, - me)] (1 - C) 



NX 



(All) 



[uj - N{kc - iKc)] [iu - lua - (N - l){kc - mc)] ' 
Similar to Eq. ( ^.14| ), the propagator can be written as a sum of two lorentzians, 



(A12) 



[uj - N{kc - i/tj] [uj - loa - {N - l){k^ - m^)] - A^A 

1 - ) 1 - aL^) 



UJ — fiV^'' UJ — 



(TV) 



(A13) 



where 



1 ± 



(^a - fee + mc)/2 

^{uJA -k, + iKc)y4: + NX 



(A14) 



(ujA-h 



C I ^'^c 



+ A^A, 



(A15) 



By substituting A^ = 0, we have A\' = 0, A[ 
the propagator reduces to 



1, fi^"* = uja — kc + iKc, fill'' = 0, and 



(0) 



(A16) 
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The propagator 



cu-H 



|i;(^-i)c) 



can be evaluated by the usual commutation relations, and noticing that 

1 1 



V 



uj-Hq ^(jv - 1 
1 



4)^-V+|O;0) 



N 



It can be proved that 



n / dhg*ih)al\0; 





u - ua - {N - l){kc - iKc) 

1 



1 



(iVO) 



(A17) 



(A18) 



(A19) 



The calculation of 



$ffH = (l;(Ar-lX 



|0;A^c: 



is similar to that of ^Q^\oi)). In fact, the symmetries in the commutation relations result in 
the relation 



(A20) 



(A21) 



The above equation can also be obtained immediately by replacing g{k) with g*{k), which 
does not change the expression in this case since the quasi-mode propagators do not contain 
any free g{k) or g*{k) terms. 
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d. <J>if (-) 



The propagator 

<l>!f H = (1; {N - l)c|^-T7|l; {N - l)e) (A22) 

to — rl 

can be evaluated by similar derivations, yielding 

(1; (AT - 11; {N - l\) = TTF—^t. ' (A23) 

and for an even integer i, 

A^A 1 / V V'"^ 

Similarly all paths with odd number times of interactions have null contributions to the 
propagator. Hence, we have 

1 



NX 



X n + 



UJ - N{kc - ?Kc)] [UJ - UJa - (N - l){kc - iKc)] - NX 
uo — N{kc — inc) 



" oM + ■ ^^^^^ 

A very simple relation exists between the e (excited-state) e and (ground-state)— g 
propagators, 

UJ - UJa - [N - l){kc - iK^) 



It can be verified that (^^iP{uj) reduces to that of Eq. ( [4.14|) if we take = 1. 
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FIGURES 

FIG. 1. A sketch of the system: A two-sided Fabry-Perot cavity with a two- level atom inside 

and partially reflecting mirrors at both ends. We only consider identical mirrors and the atom 
being locating at the center. 

FIG. 2. The left-luminating modes, with lumination from the left, and transmitted waves only 
at the right. This set of modes are labeled by the subscript L and the positve wave-number 
k. Together with the right-luminating modes, these continous field modes form a complete and 
orthogonal set of the system. 

FIG. 3. Basic components of the Feynman diagrams: (a) The total propagator from initial 
state ijji to final state ijjf. (b) Free propagation of a ground-state atom, (c) Free propagation of an 
excited-state atom, (d) Free propagation of a k-th mode photon, (e) An excited-state atom decays 
into ground state and emits a fc-th mode photon, (f) A ground-state atom excited by a k-th mode 
photon and jumps to the excited state. 

FIG. 4. The propagator $ff(a;) = g'i^ {iv; cf), cf)) can be expressed as the sum of an infinite 
series, corresponding to all the possible paths. Bold-wavy lines are used to represent intermediate 
quasi-mode photon state. This series is exactly summable. 

FIG. 5. Relationship between the propagators of single-excitation. The G's appear in different 
components are different because they are attached with different input and output legs. 

FIG. 6. Feynman diagrams for gS? {uj;ki,k'i). The first dia gram is the unlinked diagram with 
a spectator photon. The second diagram corresponds to S = {ki, k[}, which means the atom emits 
the output photon prior to the absorption of the input. All other diagrams belong to the third 
group, corresponding to S = {k[,ki}. 

FIG. 7. Relationship between the propagators with two excitations N = 2. 
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FIG. 8. Contour-plot of \C{ki, k2)\'^ for A = O.Ik^, with four different widths of the input 
photon: (a) Km = 107sp, (b) Kin = 7sp) (c) Kin = 0.57sp and (d) Kin = 0.l7sp. The axis labels are 
in unit of Kc- 

FIG. 9. The contributions of different Feynman diagrams for the case Kin = 7sp in Fig. |8| (b): 
(a) The unlinked diagram, (b) The lowest order linked diagram plus the unlinked diagram, (c) All 
diagrams. The axis labels are in unit of Kq. 
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